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Abstract 

Radiation from acoustic sources located inside randomly-layered structures is studied 
using the transfer matrix method. It is shown that in contrast to the periodically 
layered cases where the radiation can be either enhanced or inhibited depending on the 
frequency, and the characteristics and the material composition of the structures, in the 
random structures the radiation is always inhibited. The degree of inhibition depends 
on the acoustic frequency, number of random layers, and the randomness and acoustic 
parameters of the structures. Both spherically and cylindrically random structures are 
considered. The results point to the possibility of designing sonic waveguide devices 
which will not suffer from the energy loss caused by radiation, thus allowing effective 
energy confinement or long-range energy propagation. 

PACS number: 43.20.Fn. 



INTRODUCTION 

When placed in spatially structured media, the radiation or transmission of optical or acous- 
tic sources will be modulated, a fact of both fundamental importance and practical signifi- 
cance. The structure-modulated transmission, called waveguide propagation, is the backbone 
of modern opto-electronics and acousto-optics systems. Designing proper waveguide devices 
that can convey information without or with little energy loss thus has been and continues to 
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be a prime motivation for theoretical studies of wave radiation and propagation in spatially 
structured media [I]]. 

Much effort has been focused on effects of metal and dielectric interfaces, which can be 
constructed in either planar, cylindrical, or spherical forms, on optical transmission and radi- 
ation 0, |3] ^, [J. The optical transmission or radiation in periodic structures has attracted 
particular attention in different areas of applied physics, as in the periodic situations the in- 
teraction between propagating waves and structurally periodicity can be either constructive 
or destructive, leading to significant enhancement or inhibition respectively. In these situ- 
ations, the waveguides act as a filter which selects particular frequencies for propagation. 
Understanding of optical propagation in periodicity has been vital to the design of optical 
devices including optical fibers, semiconductor lasers, hetero-junction bipolar transistors, 
quantum well lasers, filters, and resonators fl], ||, [7], [|, |9|, [K| [11] . 



Propagation of acoustic waves in spatially structured media has also drawn attention 
recently Jl2|, [13], [L4], [15| . The investigation of acoustic counterparts not only paves the 
way for the possible design of new acoustic devices, but the acoustic models themselves are 
advantageous in a number of situations. First, as they are of scalar nature, the acoustic 
waves are relatively simple to handle yet not to compromise the generality, making them an 
ideal system for understanding more complicated situations with vector wave propagation. 
For instance, the recent results on acoustic propagation in water with parallelly placed air- 
cylinders make it possible to study the ubiquitous phenomenon of wave localization |16| 



in an unprecedented detailed and manageable manner [15, 17]; the localization of optical 



waves has posed a long standing problem and a subject of much debate [0. Second, the 
research is motivated by potential applications in acousto-optic fiber devices fTP[|. Third, it is 
relatively easy to manufacture hetero-structures with large contrast in acoustic impedance. 
This allows not only the study of strong scattering, but the usage of the properties of the 
strong scattering in such situations as noise reduction 



In the previous paper , the acoustic radiation from a source located inside periodically 
layered spherical cavities has been considered. It is found that significant enhancement or 
inhibition on the radiation is possible by varying the acoustic parameters and the periodicity 
of the structures of the guides. The analysis predicts well-defined peaks and nodes in the 
cavities. The fact that wave transmission in the periodical situations is only possible for 
certain ranges of frequency is useful as it is of help for devising apparatuses such as filters. 
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For such applications such as energy transport, however, it is desired that no or little energy 
be radiated for any frequencies. In other words, devices are designed so that energy are 
confined inside the devices. In this article, we consider acoustic radiation from acoustic 
sources located inside randomly layered structures. The structures can be of spherical and 
cylindrical geometries. We show that the radiation is inhibited for all frequencies for any 
given randomness, and waves are confined in the structures. The confinement or localization 
effect becomes increasingly significant as the randomness increases or the acoustic contrast 
increases for a given randomness. These features hold for both spherical or cylindrical 
structures. The results may be of help in the design of wave processing devices such as 
resonators or efficient energy guides. 

I. FORMULATION OF THE PROBLEM 

We consider a unit acoustic source located insider a layered structure, which can be either 
spherical or cylindrical geometry. The conceptual layout is sketched in Fig. 1. The most 
inner and outer radii of the structure is R\ and Rn respectively. Between Ri and Rn, there 
are N — 1 randomly placed interfaces. The boundary at i?i (i — 1, 2, . . . , N) is denoted as 
the ith interface. These N — 1 interfaces are randomly located between R\ and R N . The 
further parameters are as follows. The sound speed and the mass density inside Ri are c\ 
and pi, while the sound speed and mass density of the surround medium, i. e. outside Rn, 
are denoted by c and and p respectively. The sound speed and mass density between the ith 
and the (i + l)th interfaces are c^+i, pi + ±. We define gi = pij p and hi = q/c. The parameters 
g,i and Cj are called acoustic contrast parameters. 

The Helmholtz wave equation inside the structure is 

(V 2 + ki)p(f) = -47r<5 (c °(f), for r < R u (1) 

where V 2 is the usual Laplacian operator, k\ is the wave number {k\ = iv/ci), and 5 is the 
Dirac delta function representing the source. The superscript d refers to the dimension. For 
spherical structures, d — 3, while for cylindrical structures d = 2. The wave equation outside 
the structure is 

(V 2 + k 2 )p{r) = 0, for r > R N , (2) 
with k — uj/c. Similarly, the wave equations for the random media between R\ and Rn are 

(V 2 + k^pif) = 0, for i = 1, 2, 3, . . . , N - 1, and Ri < r < R i+1 , (3) 
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with ki = uj/c{. 

Before going into details of the problem we make a detour to consider a more general 
transmission through an arbitrary interface. The general solution to wave equation such as 
those in Eqs. ([!]), (0), and ([3]) can be written as 

p i (r)=A i G^(k i r) + B i G^\k i r), (i = 1,2,...), (4) 

where "f" means taking the complex conjugate, Ai and Bi are coefficients to be determined 
by boundary conditions, and G^ is the Green's function in the d-th dimension and is written 

as 

mH^ihr), d = 2 
d = 3. 



G^ d \kii 



1 kiT ' 



(5) 



Note that the Green's function G^ also represents the wave transmitted from the unit source 
without the presence of the layered structures; its complex conjugate represents the inward 
moving wave. 

To solve for the unknown coefficients Ai and Bi, we invoke the boundary conditions that 
require the pressure field and the radial displacement be continuous across the interfaces. 
Consider an arbitrary interface at R. The wave on the inner (<) and outer (>) sides of the 
interface are respectively denoted by 

p <t> (r) = + fl^G^r). (6) 

The coefficients A and B on the two sides of the interface can be related by the usual 
boundary conditions, 

A < G^{k < R) + B^^i^R) = AyG^ikyR) + ByG^^R), 



and 



P< 



A^'i^R) + B < G {d)] \k < R) 



-2l 



where "/" means the derivative, e. g. G^ (x) 
in the matrix form 



d_ 

At 



B. 



T(R) 



(7) 



A > G^{k > R)' + ByG^'ikyR) 

G^ d \x). Eqs. (0) and (g) can be written 

(9) 




with 
T{R) 



G^{k<R) G {d) \k < R) 



G^ikyR) G^\k>R) 



G^\k < R) ^G^'i^R) 



G^'ikyR) ^G^'' {k>R) 



(10) 



The matrix T is called the transfer matrix. 

For a system consisting of multiple interfaces, the transmission and reflection coefficients 
can be related through a consecutive product of the transfer matrices at all the interfaces. 
We denote the resulting matrix by M. Therefore, the waves inside and outside the layered 
structure are related through 

/ 



Br 



M 



B ou t 



'ID 



where 



N 



M=nw 



i=l 



mu m 12 



y m 2 i m 22 



Now we come back to the problem. Inside the most inner interface, the transmitted wave 
from the source is subject to reflection from the interface at R\. The total wave can be 
expressed as 

p( r ) = G^ihr) +p R {r), for < r < R u (12) 

where the first term is the transmitted wave from the source, and the second term is the 
reflected wave. Since the reflected wave must be finite at the origin, it can be written as 



Pr(t) = -QG^(hr) + QG^\hr). 



(13) 



The negative sign in front of the first term asserts that the reflected wave remains finite at 
the origin. Thus the total wave inside the layered structure is 



p(r) = (1 - Q)G {d) (hr) + QG {d) \k L r), 



(14) 



which gives Ai n = 1 — Q and B{ n = Q. Another observation is that there is no reflected wave 
outside the layered structure, i. e. beyond R^, we have 



B 



0, for r > Rm- 



out U ) 

Taking these into consideration, Eq. (|TT|) becomes 

1 - Q | _ ( m n m 12 
Q J y m 2 i m 22 

This equation yields the solutions 



(15) 




(16) 



^From these solutions, the radiated acoustic intensity can be computed from 

loutir) 



1 /2n\ 3 - d \A n 12 



1 out 



2pc V k ) r d ~ l 

The reflected intensity is 



W = ^-m Sr- (19) 



]_ /27r\ 3 - d |_Q| 2 
2p 1 c 1 \kiJ r 

^From Eq. (|14]) states that the effective source transmission intensity is 



1 /2n\ 3 - d \1-Q\ 2 

7 -w = 2^(17) h^- (20) 



Energy conservation states 



i-QI 2 = f-r) -MUI 2 + (1-) — \Q\ 2 - (21) 



7 / n> I 7 / (JUL 1 I 7 i 

«ol / P1C1 \ k / pc \Ki J P\C\ 



We note some typographical errors about Eqs. (21), (22) and (23) in fll4 |. 



We define the transmission (TR) and reflected (RF) coefficients as follows 

TR=|A ou< |, RF = (22) 

It will become clear that when randomness is added, the transmission will be inhibited for 
all frequencies. The energy will be confined inside the layered structure. With the random 
layers, the energy flow in the radial direction decreases as the randomness or the number of 
layers increases, a useful property for random layered cavities to transport energies. 

II. NUMERICAL RESULTS 

The situation that the coating layers are periodically placed for spherical cavities has been 
considered in [I^j. For the reader's convenience, we re-plot the results in Fig. |2|. In the 



computation, the layered structure is constructed as A — M — W — M — W — ■ • • — W, 
where A represents the air inside the cavity, i. e. the air fills the space r < Ri, M refers 
to the coating material whose acoustic impedance is p m c m , and W refers to another coating 
material and we assume it to be similar to that of the surrounding medium, i. e. the water. 
We define g = gi — p m j p and h = h 2 = c m /c. Here we take the acoustic contrasts as 
Pm/p = 2, c m /c = 2. The parameters for air are pi/p = 0.00129, C\/c = 0.23. The horizontal 
bars represent the interfaces. There are 30 interfaces. The thickness of each layer is set to 
be identical and equals 0.2i?i; thus the total thickness of the layered materials is 29 x 0.2i?i. 
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Figure ^ indicates the following. (1) Periodically layered structures selects particular 
frequencies for transmission, i. e. at these frequencies, the transmission is greatly enhanced. 
(2) The spectral valleys in which transmission is greatly inhibited are equivalent to the for- 
bidden bands observed in regular lattice solids. These valleys are called forbidden frequency 
bands. It is known from the previous results that these forbidden bands are caused by Bragg 



reflection for a periodic structure |14|]. (3) The reflection is also significantly enhanced at 
certain frequencies, and the separation of the reflection peaks is almost constant. (4) It is 
interesting to see that the results for both spherical and cylindrical geometries are qualita- 
tively similar. The transmission and reflection peaks are only shifted slightly. In the later 
discussion, we will focus on the results for the cylindrical geometry while the relevant results 
for the spherical structure will be added only when needed. 

Now we consider the randomly layered cases. In the simulation, the layered structures are 
constructed in the similar way as in the periodical structures, except that the interfaces are 
randomly placed. To be explicit, the layered structure is A — M — W — M — W — ■ • • — W. 
There are N interfaces, i. e. there are N horizontal bars. When there is no randomness, the 
interfaces are equally spaced along the radial direction; the distance between the interfaces 
is D = 0.2R\. In this way, the z-th interface is located at r« = Ri + (i — 1)D with i ranging 
from 1 to N. The level of randomness is controlled by the degree of allowing the interfaces to 
shift from their locations for underlying periodically structures. We define the randomness 
as A = ^, where S is the range within which the interfaces are allowed to shift from their 
locations when there are no disorders. For example, the location of the z-th interface can be 
randomly varied with the range between R\ + (i — 1)D — 5 and Ri + (i — 1)D + 5. Clearly, 
the total random is the case that A = 1. 

The effects of the level of randomness, numbers of random layers, and acoustic contrast 
on the radial transmission are shown in Fig. ^. As both spherical and cylindrical geometries 
have the similar features, we only plot the results for the cylindrical structures here. The 
important message from the figure is that when the randomness is added, the transmission 
is repressed for all frequencies. In other words, no energy can escaped and all energies 
are confined inside the structure. In particular, Fig. ||| shows the following. (1) When the 
randomness is introduced, the transmission becomes subdued; when the randomness is small, 
the band effects from the underlying periodic structures are still noticeable for low frequency 
bands. This is shown by the case of A = 0.3. (2) The inhibition is more significant for high 
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frequencies. For the low frequencies, the inhibition seems still possible. This is due to the 
fmiteness in the number of layers. With increasing number of layers, the inhibition will be 
extended to low frequencies. (3) Fixing the number of layers and when the randomness is 
greater than a certain value, the effect of the variation in randomness becomes not prominent 
for high frequencies. This is shown by the tendency that the curves for A = 0.3 and 1 merge 
at high frequencies, referring to the case in, for example, (c) and (e). (4) Increasing the 
number of random layers, the transmission will be reduced further. In fact, the transmission 
will decay exponentially with the number of the random layers, as will be shown later. (5) 
The above features hold when the acoustic parameters vary. (6) For the spherical geometry, 
we obtain the similar results as in Fig. ^ and because of this we do not show the results 
here. 

With increasing the number of the random layers, the transmission will decrease expo- 
nentially. This is illustrated by Fig. |] for the case g = h = 2. Here the results are averaged 
over 200 random configurations. Here is shown the transmission versus the number of ran- 
dom layers for four frequencies and four randomness levels. Out of the four frequencies, two 
are located at where the transmission is possible when no randomness is introduced, and 
one is within and one is within but close to the edge of the forbidden band of the under- 
lying periodic structure. It is clear that the increasing randomness decreases gradually the 
transmission in the regimes in which the transmission is possible when there is no disorder. 
These regimes may be called the passing bands for the underlying periodical structures. This 
property is also true for frequencies with but in the vicinity of the forbidden band. Inside 
the forbidden bands, however, though inhibited the transmission is increased with the added 
disorder as compared to the case without disorders. This property has also been observed 
in one-dimensional random liquid media [f2~T|, [22]] . This indicates that the inhibition due to 
randomness and the inhibition due to Bragg reflection may have different causes, and the 
two effects compete. For all frequencies, be within the passing or forbidden bands, the trans- 
mission decreases exponentially with increasing numbers of random layers, implying that the 
most energy is localized near the transmitting source. 

The relation between the transmission and the number of layers can be approximately 
described as 

' = ^- m ' N '- (23) 
The parameter N represents the effective number of random interfaces to localize the energy, 
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and is named the localization length in terms of the number of interfaces. The localization 
lengths for the cases in Fig. ^ are summarized in Table [I]. It is clear from the table that within 
the passing bands of the underlying period structure the localization length decreases with 
increasing frequency and disorders, while well within the forbidden bands the localization 
length decreases with increasing disorders. This table also shows that the localization lengths 
are almost identical for the spherical and cylindrical geometries. Moreover, we also see that 
for the totally random configurations, the energy confinements can be achieved by just a few 
random interfaces. Further simulation indicates that the localization length decreases with 
increasing the acoustic contrast. 

Table 1: Localization length vs kR\ for various randomness A. 

Spherical Shell Cylindrical Shell 

kRi\A 0.1 0.3 0.5 1.0 0.1 0.3 0.5 1.0 

5.5 355.7 42.9 18.0 7.7 390.3 44.2 17.9 7.7 

6.0 145.2 20.8 11.0 6.2 143.7 20.4 11.1 6.3 

6.5 12.4 8.8 6.9 5.5 12.4 8.8 7.0 5.4 

8.0 2.0 2.1 2.7 4.5 2.0 2.1 2.7 4.5 



We also examine the reflection behavior. The conservation law in Eq. (|2l|) is rewritten 

as 

When the radiation is stopped, i. e. \A out \ — > 0, we are thus led to the relation 

|i-Q| 2 ^|Q| 2 . (25) 

The numerical confirmation is shown in Fig. [5], in which the ratio |1 — <5|/|Q| is plotted as 
a function of frequency in terms of kR\. With reference to Fig. |2|, we see that the ratio 
deviates from one only for frequencies within the passing bands for the regular structures. 
With added disorders, this ratio is virtually one for all frequencies. Note that the reason why 
the ratio is close to one even for the periodic structures is that we take out a term associated 
with the acoustic impedance. 

Finally we note that the results in this papers bear some similarities to the Anderson 
localization in one dimensional random systems in that no waves can propagate in such a 



systein [ ZI I . The most significant difference is that in the present case, the waves are localized 
near the transmitting site. In the ID cases, however, the energy needs not be confined near 
the source, and there is a stochastic resonance behavior which is absent from the present 
situations. 



III. SUMMARY 

In this paper, we consider acoustic radiation from randomly layered cavities. The results in 
the present paper convey the information that when a cavity is coated by random layers, 
virtually no energy can be radiated in the radial direction. The waves are mainly localized 
inside the cavity, and its decay along the radial direction follows an exponential law. The 
results presented here may be useful for designing acoustic 'lasers', resonators, or energy 
transport. For example, suppose there is a cylindrical wave guide. The propagating wave 
can be generically written as e ikzZ f(r) in the cylindrical coordinates. When the guide is 
coated by random materials, according to the results, no energy can be radiated into the 
transverse directions, all energies will be confined inside the guide and can only propagate 
in the longitudinal direction, i. e. f(r) decreases with r. By adjusting the coating materials, 
the energy transport along the longitudinal axis may be tuned to fit the applications. The 
results may also be useful for designing possible 'sonic fibers' in analogy with the recent 
all-dielectric optical fibers |J. For the spherical waveguides, various energies can be stored 
inside the guides by adjusting the coating materials. 
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Figure 1: Conceptual layout for spherical and cylindrical cavities. 
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Figure 2: Transmission and reflection coefficients as a function of frequency in terms of kR\ 
for periodically layered cavities: (a) Transmission; (b) Reflection. 
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Figure 3: Transmission as a function of frequency in terms of kR\ for randomly layered cylin- 
drical structures for various disorders, acoustic contrasts, and numbers of random layers. The 
results shown are for an arbitrary random configuration. 
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Figure 4: Transmission as a function of the number of random layers for various frequencies 
and different randomness 
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Figure 5: The ratio |1 — Q\/\Q\ versus kR±. The solid and dashed lines are for the periodic 
and random structures respectively. Here g = h = 2 and iV = 100 
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